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On the Diagonalization of 
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J. F. van Diejen 



ABSTRACT. We discuss the simultaneous diagonalization of a family of com- 
muting difference operators by Koornwinder's multivariable generalization of 
the Askey-Wilson polynomials. The operators constitute a complete set of 
quantum integrals for a difference version of the n-particle quantum Calogero- 
Sutherland system related to the root system BC n . 



1. Introduction 

Some years ago, Ruijsenaars introduced [19] an integrable quantum n-particle 
model that can be viewed as a relativistic generalization of the well-known quantum 
Calogero-Sutherland system [22, 23]. The commuting quantum integrals for the 
model are given by difference operators of the form 

(1.1) 

H r = ~[iv*(xj - x k ) e (3 ^j£ j9j ~[[vi(x k - xj), r = l,...,n, 

Jc{i,...,n} jeJ jeJ 

I J I izz t k ^ J Jz ^ J 



where dj = t£-, | J| denotes the number of elements of J C {1, . . . , n}, 
(1.2) v(z) 



sin f (i/3g + z) 



sin(^2;) 

(One has e tl3d ' ^(xi, ... , x n ) = ^(xi, . . . , Xj-i , Xj + if3, Xj +1 , ... , x n ), so H r is a 
difference operator of order r.) 

For g = the coefficients of H r (1.1), (1-2) are trivial (v = 1) and the difference 
operators are simultaneously diagonalized by the (Fourier) basis of symmetrized 
exponentials (or monomial symmetric functions) 

(1-3) 

m x (x) = e ia £"=i A ^ 3 ', A e A = {A <E Z™ | Ai > A 2 > • • • > A„ }. 

XES„\ 

For arbitrary g > a basis of joint eigenfunctions of Hi, . . . ,H n follows from 
work by Macdonald on orthogonal g-polynomials in many variables [12, 13, 14]. 
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One has (assuming a, (3 > 0) 

(1.4) H r * A = E r (e- a ^ Xl+f "\... , e -«/5( A "+^)) * A) A £ A, 
with 

(1.5) E r (t u ... ,t n ) = II *j » Pi = |s("+l-2j), 



JC{1,... ,n} jGJ 
J|=r 



and eigenfunctions of the form 
(1.6) * A = A 1 / 2 p\, A= n 



(„-a/3g„ia(xj-x h ). „-a/3\ ' 

where pa denotes a trigonometric polynomial that is determined by the conditions 

i. px = m A + £ c a ,a' m A ', c a ,a' £ C; 

A'GA,A'<A 

ii. (pa,™a')a=0 if A' < A. 

Here {•, -)a is the L 2 inner product with weight function A (1.6), where, as usual, 

oo 

(a;g)oo = l[(l-aqr). 

p=0 

The (partial) ordering of the cone A (1.3) is defined by A' < A iff ^'j ^ S^Li -\? 

for m = 1, . . . , n, with equality holding when m = n (A' < A iff A' < A and A' ^ A). 
Clearly, p\ amounts to m\ minus its orthogonal projection with respect to {•, -)a 
on span{m A '}A'GA,A'<A- 

The nonrelativistic limit corresponds to sending the difference step size (3 to 
zero. Apart from a constant factor the weight function A (1.6) converges to 

(1.7) A = JJ I sin f (a, - x k )f 9 , 

l<j<k<n 

1/2 

and the eigenfunctions "J" a (1-6) go over in "J'a.o = A p\,o, where £>a,o denotes the 
trigonometric polynomial determined by conditions L, ii. with the weight function 
A (1.6) replaced by A (1.7). In order to study the transition (3 — » for the 
commuting quantum integrals it is convenient to introduce [20] 

(1.8) H r = (3~ r £ r = l,...,n 

0<p<r ^ T ' 

(taking Ho = 1 by convention). The operators H±, . . . ,H n generate the same alge- 
bra of commuting difference operators as H±, . . . , H n , but they have the advantage 
of having a nontrivial limit when the step size goes to zero. From Eqs. (1.4), (1.8) 
it follows that 

(1.9) 

H r * A = P~ r E r (1 - e - al3{M+ ^\. . . , 1 - e -^( A -+^)) * A , A € A, 
which leads for (3 — * to 

(1.10) H rfi * A ,o = a r E r (Ai + p u . . . , A„ + p n ) * A ,o, A <= A, 
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with 

(1.11) # r , = lim£T r = (-i) r IJdj + l.o. 

JC{1,... ,n} j£J 
\J\=r 

(where l.o. stands for terms of lower order in the partials dj). Explicit computation 
of H rt o for r = 1, 2 yields 

(1.12) 

l<j<n l<j<fc<ra 2 V J ' 

The differential operators H 10 ,...,H nQ constitute a complete set of commuting 
quantum integrals for the nonrelativistic Calogero-Sutherland system. Explicit ex- 
pressions for the integrals can be found in Refs. [3, 21] and [16]. 

It is well-known that the Calogero-Sutherland model admits a Lie-theoretic 
generalization in terms of root systems [17]. From this viewpoint the models in 
this introduction correspond to the root system A n -\. In the rest of the paper 
we are going to focus on the more general case corresponding to the root system 
BC n . Specifically, we will discuss the diagonalization of a complete set of quantum 
integrals for a difference version of the quantum Calogero-Sutherland model related 
to the root system BC n . This difference quantum model may be regarded as a 
BC„-type generalization of Ruijsenaars' relativistic Calogero-Sutherland system. 

Remarks: i. The polynomials p\ were introduced by Macdonald [12]. They 
form a g-version of the so-called Jack polynomials [11] (with q = exp(— af3)). The 
transition p\ — » p\ t o (/3 — > 0) corresponds to the limit q — » 1. In this limit Mac- 
donald's polynomials reduce to Jack polynomials. 

ii. Conjugation of H r with A 1 / 2 eliminates the square roots in the coefficients 
of the difference operators: 

(1.13) A -1 / 2 H r A 1 / 2 = nv(x j -x k )e if> 'Zie' Bi . 

JC{1,... ,n} jeJ 

\J I — v k^J 

The transformed operators (1.13) are diagonalized by the Macdonald polynomials. 
They are essentially the same as the Macdonald difference operators associated with 
the fundamental weights of the root system A n -\ [13, 14]. (The exact connection 
is established after restriction to the center of mass plane x\ + ■ ■ ■ + x n = 0.) 

2. Commuting difference operators 

The quantum integrals for our model are given by 

(2.1) H r = Uj^jiVtfj.eW-'V^j.j., r = l,...,n, 

JC{1,... ,n}, | J|<r 

Ej=±i, je J 
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with d eJ = Ej G J £ i 5 i' 



(2.2) V sJ . K = n w(£jXj) JJ v(£ j x j +e j .x j ,)v(e j x j +e j .x j , +ip) 

3^ J 3,3 'G J 

j<j' 



X H "(^j + X k )v(£jXj - X k ), 



jeJ 

keK 



(2-3) 

^if,m = (-l) m ^ Qw(£;2!;) j| v(eixi + £ V X V )v(-£lXl - £ V X V -ij3) 

LcK, \L\=m leL 1,1'eL 
£i = ±l, *Gi l<l' 

x j| v(eixi + x k )v(sixi - x k ), 



lEL 
k£K\L 



and 

u(z) 



sin ^(iftg + z) 



sin(|-z) 
(2.4) 



sin f (ipgg + z) cos f + z) sin f + f + z) cos |(jgg| + f + z) 
sin(fz) cos(fz) sin |(i£ + z ) cos |(^ +z ) 

It is clear that, just as the quantum integrals for the models considered in the in- 
troduction, these difference operators are invariant with respect to permutations 
of the particle positions x\,... ,x n . However, the operators are no longer trans- 
lational invariant; instead, we now have invariance with respect to the reflections 
Xj — » —xj, j = 1, . . . ,n. Furthermore, apart from the potential v, which governs 
the interaction between the particles, there is now a second potential w modeling 
an external field. From now on we will always assume that 

(2.5) a,p>0 and ff,^ ,^ > 0. 

These constraints ensure that our quantum integrals H r (2.1)-(2.4) are self-adjoint. 
(One has v(xj) = v(-Xj), w(xj) = w(-xj), so V eJ - K = V^ eJ - K and U K , m = 

UK,m-) 

The (weight) function that transforms away the square roots in the coefficients 
of the difference operators (cf. Remark ii. in the introduction) now reads 

(2.6) A = Q d v (ex,j + c' xjf) jj[ d w (sxj), 

1<3<3 <n l<j<n 
e,e' = ±1 e = ±1 



with 
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(e iaz ;e- a ^) 00 (-e iaz ; e~ a ' 3 ) oc , 



d w (z) 



' (e- a ^o e iaz. e-a/?)^ (_ e -«/3gi giaz . g-a/3 ) ^ 

( e ia(i/3/2 + z). (_gia(i/3/2+z) . e -<*0) c 



( e -a/3g' a gia(i/3/2 + z). (_g-a/3fli gia(i/3/2 + z) . g-a/3 ) ^ ' 

Conjugating with A 1 / 2 results in 

(2.7) 7i = A-^F.A 1 / 2 = 5] ^,,-|j|^ J;J c e ^-. 

JC{1,... ,n}, J|<r 
£j=±l, iG J 

The commutativity of Hi,... ,H n (or equivalently 7i\,... ,1~L n ) is by no means 
immediate from the above expressions, but it follows from the fact that (as we will 
see below) the operators have a basis of joint eigenfunctions. 

3. Triangularity 

The first step in the diagonalization process consists of demonstrating that the 
transformed difference operators T~L r (2.7), (2.2)-(2.4) map the space of permutation 
invariant, even, trigonometric polynomials into itself. This space is spanned by the 
basis of W-invariant monomials 

(3-1) 

mx (x) = e ia ^=i X ' jXj , A G A = {A G Z" | Ai > A 2 > • • • > A„ > }, 

A'GWA 

where W is the group generated by permutations and sign flips (W — S n x (Z2)"). 
When acting with Ji r on a monomial m\ we obtain a W-invariant trigonometric 
function that is rational in exp(±iaxj), j = 1, . . . ,n: 

(3.2) ii r m\ = ^ V eJ+iE j_ ; jc nJ c m A (x + i/3 J2 j€J+ £ j e j) > 

J + nJ_=0, J + UJ_ \=r 

with ej denoting the jth unit vector of the standard basis of ffi n and 
(3.3) 

V e j +t£ j_ ; K = (-1)' J_ {[ w( £ j x j) ]| v(ejXj +x k )v(e j x j - x k ) 

jGJ+UJ_ jGJ+UJ_ 
kGK 

X JJ ^(SjXj + Ej-Xj, ) ^(SjXj + Ej-Xj, + if}) 
3<3 

x Y[ v ( e j' x j + -•./'•'•./' I ''I \i- r .i ~ £ j' x i' ~ ifi) 
jJ'eJ- 

3<3 

X Y[ V ( £ 3 X 3 + X j') V ( £ 3 X 3 ~ X 3')- 
3&J+ 

j'eJ- 

To prove that 7i r m\ is actually polynomial in exp(±iaa;i), . . . , exp(±iax rl ), it 
suffices to demonstrate that, viewed as a function of Xj (with all other variables 
fixed in general position), the r.h.s. of (3.2) is regular. Thus, we must show that 
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the (generically) simple poles in the terms of 7i r m\, which are caused by the zeros 
in the denominators of the potentials v,w (2.4), all cancel each other in (3.2). The 
poles are congruent (mod 2ir/a) to one of 

x j ~ (+ f)> type !> 

Xj ~ ±x jf , j' ^ j, type II, 

~ ±i(3/2 (+ type III, 

Xj ~ ±i/3 ± aij- , j' ^ j, type IV. 

For the poles of type I and type III it is enough to consider only those congruent to 
and ±i/3/2, respectively. This is because a translation over the half-period ir/a 
is equivalent to an interchange of parameters <-> g['^ . Let us first consider the 
types I and II. The residue at Xj = vanishes because (3.2) is even in Xj. The 
vanishing of the residues at the poles of type II is also immediate from the W- 
invariance: we first use the permutation symmetry to conclude the vanishing of the 
residue at Xj = xji , and then the invariance under sign flips to conclude the same 
for the residue at Xj = —xy. For the remaining two cases (type III and type IV) 
we may assume that all signs are negative. This corresponds to considering terms 
in (3.2) with ej, £y = +1. The poles of type III arise in the terms with j either in 
J + or in J . By combining each term with the term corresponding to the index 
sets obtained by bringing j from J + to J_ or vice versa, it is seen that the residues 
cancel each other pairwise. The situation for poles of type IV is very similar; these 
poles only occur in terms with both either in J + or in J_. By bringing j, j' 
from J_|_ to J_ or vice versa, one infers that the residues of the corresponding terms 
also cancel in pairs. 

Now that it has been verified that the r.h.s. of Eq. (3.2) is regular in Xj, j = 
1, . . . ,n, we may conclude that T~L r m\ is a B^-invariant trigonometric polynomial. 
Hence, it can be written as a finite linear combination of the monomials (3.1). It 
turns out that the expansion of TL r ra\ in monomials is of a very specific form: 

(3.4) H r m x = [tir]\,y my, with [K]>, A ^C, 

A' 6 A 
A'<A 

where the (partial) ordering of the cone A (3.1) labeling the monomial basis is 
defined by 

(3.5) A'<A iff X 'j ^ J2 X i for m = l,...,n. 

In other words: the operator T~L r is triangular with respect to the partially ordered 
basis of B^-invariant monomials {toa}aga- 

The triangularity of Ji r hinges on the fact that the potentials v(z), w(z) (2.4) 
have constant asymptotics for Im(z) — > ±oo: 

(3.6) lim v(z) = e ±Q/3fl / 2 , lim w(z) = e ±a ^ 9o+in+9 ''> +9 ^l 2 . 

Im(z) — >±oo Im(.z) — *ioo 

To see this we put x = iRy, with y € ffi™ and y\ > y 2 > ■ • ■ > y n > 0, and combine 
the asymptotics 



(3.7) my (iRy) = «""£'.' A ' " J (1 + o(l) ) for 



R 



ON THE DIAGONALIZATION OF DIFFERENCE CALOGERO-SUTHERLAND SYSTEMS 7 



with Eq. (3.6) to conclude that 

(3.8) (H r m x )(iRy) = 0(e aR ^=^ jVj ) for R -► oo. 

By comparing the asymptotics (3.7) and (3.8), and using the fact that 
(3.9) 

A' < A iff J2 K Vi ^ J2 X i y i yy e W with Vi> ■■■>Vn>0, 

l<j<n l<j<n 

one infers that in the expansion of 7i r m\ in monomials my only the terms with 
A' < A may occur with a nonzero coefficient. 

4. Eigenfunctions and eigenvalues 

Let p\ be m\ minus its orthogonal projection onto span{mA'}A'GA, A'<A with 
respect to the inner product {•, -)a (i-e., the 1? inner product with weight func- 
tion A (2.6)). Clearly p\ is uniquely determined as the monic polynomial in 
span{m-A' }a'ga, A'<A that is orthogonal to span{m-A' }a'ga, A'<a- The triangular- 
ity oiTL r ensures that TL r p\ lies also in span{m-A' }a'ga, A'<a 5 and the symmetry of 
7i r with respect to {•, -)a implies that T~L r p\ is orthogonal to span{mA'}A'eA, A'<A- 
(One has (H r px,mx')A = (px,T~Lr ?tia')a = if A' < A.) Hence, TL r p\ must be 
proportional to p\, i.e., p\ is an eigenfunction of Ji r . 

Because the basis transformation {toa}aga — * {pa}aga does not affect the 
matrix elements on the diagonal, the corresponding eigenvalues are given by (using 
Eqs. (3.4), (3.7) and (3.9), together with the fact that in the last equation the 
inequalities are strict if A' < A): 

(4.1) [7i]A,A = lim e~ aR ^"=^ iVi (H r m x )(iRy). 

R— »oo 

To evaluate this limit we use Eqs. (3.2), (3.3) and (3.6), (3.7) to obtain 



[Hr]x,x = e (-i) |j -' n ^ s ^ +pi) n 1 

j + nJ_=0, | j + uj_ \=r jeJ+ jeJ- 



a/3ej (90+Sl +g' +g[ )/2 



n 



e 



ap Sj g 



= 2 r (- 1 ) r_|J| UchafiiXj + Pj ) cha/3p h ---cha/3p lr _ lJl , 

JC{l,...,n} jeJ r<h<-<l r _ l j l <n 

0<\J\<r 

with pj = (n—j) g + (go + 9i + 9o + 9i)/2- Transforming back to Lebesgue measure 
leads to the eigenfunctions of Hi, . . . , H n (2.1)-(2.4). 
Let us summarize the final result: 

H r V x = 

E r (cha/3(Ai + pi),... ,chaf3(\ n + p„); ch (a/3p r ), . . . ,ch(a/3p n )) * A , 
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with 

E r (*!,... ,t n ; Pr ,... ,p n ) = Y (-l) r_|J| II*J Yl Ph---Ph- U \ 

JC{l,...,ra} j£J r<h<-<l r _\j\<n 

0<\J\<r 

and eigenfunctions of the form 

tf A = A 1 / 2 pa, A e A = {A e Z n I A! > A 2 > • • • > A„ > }, 

where A is taken from (2.6) and p\ denotes the trigonometric polynomial deter- 
mined by the conditions 

i. p\ = m\+ ^ cx,ymy, c a ,a' G C; 

A'eA,A'<A 

ii. (pa,™a')a=0 if A' < A 

with the partial ordering of the cone A as defined in (3.5). 

Remark: For r = 1 the operator Ji r (2.7), (2.2)-(2.4) becomes 

(4.2) n ± = ^ w ( £X j) Y[ v ( £X j +x k )v(sxj -x k ) (e if3ed i - 1). 

l<j<n kjij 
£ = ±1 

This operator and its eigenfunctions jj^, A 6 A, were introduced by Koornwinder 
[10] as a generalization of Macdonald's construction [13] for the root system BC n . 
All Macdonald polynomials associated with (admissible pairs of) the classical root 
systems (i.e., A n _i, B n , C n , D n , and BC n ) can be seen as special cases of the 
Koornwinder polynomials {pa}aga [6]. In the case of one variable (n = 1) Koorn- 
winder's polynomials reduce to Askey- Wilson polynomials [1]. 

5. Concluding remarks 

5.1. Trigonometric identities. In Ref. [6] we studied the diagonalization 
of difference operators 7i r (2.7) with V e j ; k as in (2.2) and Uk,™ given by 

(5.1) U K , m = 

( _i ) p v eLl . K \ Ll v e ( L2 \ Ll ). K \ Ll ■■ ■ v e ( Lp \ Lpl) . K \ Lp 

QCLiC-'-CLpCK, l<p<m 

(where v and w are again taken from (2.4)). We found that these difference opera- 
tors are simultaneously diagonalized by Koornwinder's multivariable Askey- Wilson 
polynomials, with the same eigenvalues as the operators in the present paper (i.e., 
with UK, m given by (2.3)). Because the difference operators are uniquely deter- 
mined by their eigenvalues [6], it follows that the r.h.s. of Eq. (2.3) equals the 
r.h.s. of Eq. (5.1). This equality hinges on nontrivial functional identities for the 
trigonometric function v (2.4) [7]. 

5.2. The limit /3 0. For (3 -> 0, one has [6] 

H r = f3 2r H r , +o(/3 2 "), H rfi = (-1)" J2 Ii d i +Lo -' 

JC{1 n}j€J 

|J|=r 
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with 

H r ,o (A}p A ,o) = 

a 2r E r ((\ 1 +p 1 ) 2 ,... ,(X n +p n ) 2 ;p 2 r ,... ,p 2 n ) Ag p x , , A € A, 

where p\ denotes a multivariable generalization of the Jacobi polynomials [4, 9] 
determined by the conditions L, ii. with weight function 

A = j[ | sin f(xj + x k ) sin f (xj - x k )\ 29 j[ | sin(f Xj)^ 90 | cos(|Xj)| 291 , 

l<j'<fc<n 

and go = go + g' , g\ = g\ + g[. The differential operators -Hi,o, • • • , H n ^ constitute 
a complete set of quantum integrals for the Calogero-Sutherland system related to 
the root system BC n . For r = 1 the operator reads 

K ° = ~Ti d > + -2^-V a2 E Ln 2 f(l J+ x k ) + S in 2 f( 1 x J -x k) ) 

l<j<n l<j<k<n 2^-7 ' 2 V 3 ' 

1 q2 y M-l) + 91& - 1 ) ) + cons tant. 

Explicit expressions for the higher-order quantum integrals can be found in Refs. 
[4, 16]. 

5.3. Differential- and difference-reflection operators. For arbitrary root 
system the existence of a complete set of quantum integrals for the Calogero- 
Sutherland system can be proved using a construction involving certain differential- 
reflection operators known as 'Dunkl operators' [9]. This approach does not yield 
explicit expressions for the quantum integrals, but it does provide so-called shift 
operators. With these shift operators well-known conjectures due to Macdonald 
for the values of normalization constants (p\,o,P\,o)a can De proved [18]. A par- 
tial generalization of this method to the Macdonald polynomials related to reduced 
root systems was presented recently by Cherednik [2]. Specifically, it follows from 
Cherednik's work that for special values of the parameters g^\ g^ the existence 
of a complete set of quantum integrals for /Hi (4.2) can also be shown with aid 
of difference-reflection operators, and that, in these special cases, Macdonald's q- 
version of the normalization conjectures for (p\,P\)a can be proved. 

5.4. Quantum group interpretations. It is expected that for certain spe- 
cial values of the parameters the Koornwinder-Askey- Wilson polynomials can be 
seen as zonal spherical functions on quantum versions of homogeneous spaces of 
classical simple Lie groups. Our commuting difference operators should then cor- 
respond to radial reductions of the Casimir elements. For the root system A n -\ 
such a quantum group interpretation of the Macdonald polynomials was given by 
Noumi for g = 1/2 and g = 2 [15]. A different approach, valid for arbitrary values 
of g, was developed in Ref. [8]. 

5.5. Elliptic generalizations. It is well-known that the trigonometric Calo- 
gero-Sutherland system has an integrable generalization involving elliptic potentials 
[17]. Explicit quantum integrals for the elliptic version can be found in [16] (differ- 
ential case, all classical root systems) and [19] (difference case, root system A n -\). 
A partial generalization to elliptic potentials of H r (2.1)-(2.4) was presented in [5]. 
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